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£H ■ Abstract 

We perform a thorough analysis of the VP-f(-f*) and ^77(7*) decays in the resonance chiral 
theory, where V stand for the vector resonances p, K* , u), <fi, P stand for 71", K, 77, rj' and 7* subsequently 
decays into lepton pairs. Upon imposing QCD short-distance constraints on resonance couplings, the 
uj — > 717(7*), P ~~ ^ K*° ~ * K°j processes only depend on one free parameter and n — > 77(7*) 

can be completely predicted. The four mixing parameters of the rj — rf system, i.e. two mixing angles 
#8,#o and two decay constants F%,Fq, are determined from radiative decays involving rj or rj '. The 
higher order low energy constants of the pseudo-Goldstone Lagrangian in the chiral anomaly sector 
are predicted by integrating out heavy resonances. We also predict the decay widths of p — > ire + e~ , 
rj —¥ 7e + e~ and </> — > rjp + p~ , which can be compared with the future measurement in these channels. 

PACS: 12.39.Fe, 13.20.Jf, 11.15.Pg 
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■ 1 Introduction 

(N ; 

t-h ■ To study the properties of rj and rj mesons is a very interesting subject in hadron physics. The reasons 
behind are twofold. First, several experimental collaborations have started or planned programs to launch 
\ the measurements of the processes involving the r\ and rj mesons with high statistics and high precision, 



such as KLOE [Ij, Jefferson Lab [2] and BES-III [3]. The huge data sample, for example 63 million events 
for r/ decays and 61 million events for rj decays expected at BES-III [3], apparently needs more and finer 
theoretical work for the analysis. Second, on the theoretical side, rj and rj mesons present important 
\ information of low energy dynamics of QCD: the mechanism of spontaneously chiral symmetry breaking 
and the Ua{X) anomaly. 

The responsible theory in the low energy region of QCD is Chiral Perturbation Theory (xPT) [1], 
whose degrees of freedom are the pseudo-Goldstone mesons, i.e. w, K, rj, resulted from the spontaneously 
chiral symmetry breaking from SU(3) L SU(3) R to SU(3)v=l+r- xPT has been proved to be a very 
successful effective field theory to describe the low energy physics of QCD [5j, that is constructed with 
respect to chiral symmetry and arranges its effective action in the expansion of momenta and the mass 
of pseudo-Goldstone bosons. Due to Ua{1) anomaly, rj is prevented to be the ninth pseudo-Goldstone 
boson. Nevertheless, when the number of the colors in QCD, Nc, becomes large, the effect of ?7a(1) 
anomaly is suppressed. Thus the rj meson becomes the ninth pseudo-Goldstone boson in the large Nc 
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limit [6] and can be incorporated into a chiral Lagrangian. Such an effective theory explicitly including the 
rf meson extends the standard SU (3) L (£) SU(3) R xPT to the U(3) L (£)U (3) R version, whose Lagrangian 
up to 0(p ) has been thoroughly investigated in [7J [8] for the even intrinsic parity sector. 

The application of the U(3) L (^) U(3) R xPT to rj — rf mixing has been performed in the literature [9, 
[TU], 111], where the mixing parameters of rj and 77' are determined by including higher order corrections, 
such as the low energy constants(LECs) and the loops. A novel finding after including the higher order 
contributions in xPT is that the conventional one-mixing-angle description for 77 — rf mixing is not valid 
any more and two-mixing-angle scheme is then proposed in [TO]. In our current discussion, instead of 
working in more detail in the top-down method to address the 77 — rf mixing problem, we are going to 
determine the mixing couplings in a bottom-up way, i.e., we assume the validity of the two- mixing-angle 
description and then directly fit them using experimental data from the relevant physical processes that 
involve rj or r/. We will focus on the radiative decay processes in the present work, since they can provide 
a large sample of data [12] that allows us to better extract the 77 — rf mixing information and they are 
less contaminated by the strong final state interaction comparing with the hadronic decay modes, such 

as rj' — > TjTTTT [13]. 

Among the radiative decay processes with 77 or rj' meson in the low energy sector, many of them consist 
of one vector resonance, such as cp — > 777, rj' — > uij, and so on [12]. Apparently, these processes are already 
beyond the validity region of xPT due to the appearance of the heavy vector resonances. It is by no means 
trivial to systematically include the heavy vector resonances in xPT, since the expansion parameters of 
%PT are no longer valid after the inclusion of the heavy multiplet of resonances. Nevertheless, the 
framework developed in Ref. [14] . named Resonance Chiral Theory (R^T), has been proven to be useful 
and may shed light on the proper construction of the Lagrangian theory that one could use to describe 
the dynamics with both pseudo-Goldstone mesons and resonances. It can be better understood within 
the framework of the large Nc QCD as theory of hadrons [15] . While in the strict large Nc QCD there 
is an infinite number of zero-width hadrons in the spectrum, in practical realization one usually needs to 
truncate the infinite tower of resonances to the lowest multiplet for each quantum number. There has 
been a large amount of research works based on this approximation, varying from determination of the 
XPT LECs pH QI] to the study of tau decays [M 03 ED] and Green functions of QCD currents [2JJ [22l 

R%T, although well respecting chiral symmetry and constructed in the guide of 1/Nc expansion, is 
still lacking of QCD dynamics at the high energy scale where the continuum is reached and perturbative 
QCD is the responsible theory. Thus it is crucial to match the behaviors of the effective field theory and 
QCD at high momentum transfer to implement as many QCD features as possible. Research along this 
line indeed has been intensively performed in many works |21[ 1221 1251 124"! 125) . This procedure directly 
results in the constraints on the resonance couplings and hence makes the theory more predictable. 

In the present work, we utilize R%T to analyze the VP^^*) and Pjj("f*) processes. Comparing with 
the work in Ref. [25], we generalize the R%T Lagrangian with the octet of pseudo-Goldstone mesons to 
the version with nonet, thus allowing us to study the processes with rj' meson. Our work is also especially 
devoted to the determination of the rj — rj' mixing parameters by fitting data, which we will explain in 
detail later in the text. 

We organize the article as follow. A mini-review on rj — rj mixing is given in Section [2J The structure 
of the relevant R%T Lagrangian is elaborated in Section [3j The computation of the decay amplitudes 
of ^^7(7*) and P77(7*) is noted in Section UJ The QCD short distance constraints are discussed in 
Section [5l Phenomenology discussion is given in Section [6] and we conclude in Section [JJ 
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2 Mini- review on 77 — rj mixing 



77 — 7/ mixing is an interesting subject in hadron physics. In the literature, the mixing angles have been 
defined with respect to different bases: the octet-singlet flavour basis and the quark flavour basis. For 
our purpose, we will always adopt the octet-singlet flavour basis to define the mixing angles throughout 
this article. The two-mixing-angle description has been proposed to settle the 77 — rf mixing |10| . going 
beyond the old one-mixing-angle description [26]. The requirement of the two mixing angles can be better 
understood in the frame. The leading order Lagrangian of U(3) L (^) U(3) R xPT is 



F 2 

+ — Ml In 2 detu, 
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The last term in Eq.([T]) represents the Ua(X) anomaly of QCD, which gives rise to the singlet 771 mass, 
Mq. F is the value of the pseudo-Goldstone decay constant in the chiral limit, with the normalization of 
Fx = 92.4 MeV. For other chiral building blocks, see Ref. [13] and references therein. 

At leading order, only one mixing angle is needed to diagonalize the octet r/g and the singlet r\\ to 
get the mass eigenstates 77 and 77', since there is only one mixing term in the mass sector. However, the 
higher order corrections, including the loop contributions and the LECs from the higher order Lagrangian, 
contribute not only to the mixing of mass but also to the mixing of kinetic term. In order to get the 
physical eigenstates of 77 and 77', one needs three steps: first to diagonalize the kinetic term, then to 
perform the normalization of each field and to diagonalize the mass term in the end, which indicates four 
parameters, i.e. two angles and two normalization constants, are needed in this procedure. It has been 
shown in Ref. |10[ [TT] the mixing can be parameterized as, 



1 

F 



F$ cos 6*g 
Fx sin Ox 



-F sin 6q 

Fq COS#q 



(4) 



where Fg and Fq correspond to the weak decay constants of the axial octet and singlet currents, respec- 
tively. By setting i*g = Fq = F and #0 = $8 i n Eq.([T]), the conventional one-mixing-angle scheme can be 
recovered. 

The works presented in Refs. |10t [TT] [9] are devoted to the determination of the mixing parameters 
Fg, Fq,9$,0q by including the higher order corrections in U(3) xPT. Recently assorted methods along 
this line have been done to settle the 77 — 77' mixing in |27[l28j. As already advertised in the Introduction, 
instead of considering more of higher order corrections to calculate the mixing parameters, we will first 
adopt the two- mixing-angle scheme described in Eq.Q for the 77 — 77' system and then determine the 
unknown mixing parameters phenomenologically. Similar works within this context have been carried out 
in Refs. [29][^[^[Sl[Ml[331[Ml[Ml[371[38l[S], which confirm the robustness of the two-mixing-angle 
description scheme. In the present work we would like to readdress the similar processes, such as the 
VPj and P77, in the framework of RxT. The advantage of R^T is that it preserves the chiral symmetry 



3 



in the low energy sector, respects the high energy behavior of QCD and incorporates all of the symmetry 
allowed operators in the construction, in contrast to the previous works where only a single constant term 
is introduced to describe the interaction vertex for VP7. In addition, we perform a global fit by including 
all of the experimental available processes with the types of V — > P7, P — > V"f and P — > 77 and also 
discuss the processes with an off-shell photon decaying into lepton pairs: V — > Pl + l~, P — > Vl + l~ and 
P — > 7/ + Z~, where V stand for p, lj, (j>, K* and P stand for n, K, rj, rf . 

3 The relevant Lagrangian of RxT 

The resonance chiral effective Lagrangian describing vector-photon-pseudoscalar and vector- vector-pseudoscalar 
vertexes with the vector resonance in the antisymmetric tensor filed formulation has been given in 
Ref. [25], in which only the octet of the pseudo-Goldstone mesons is included. Using the constraints 
derived from the short distance behavior of QCD, the Lagrangian has been used to predict the decay 
widths of oj — > 7T7 and ir — > 77, which are in good agreement with the experimental data [12] . In order 
to study the similar processes with rf , we need to generalize the existing resonance Lagrangian with the 
pseudo-Goldstone octet to the one involving the singlet state in addition to the octet. Thanks to large 
Nq QCD, this can be simply accomplished by extending the content of the unitary matrix u(x), defined 
in Eq.([2]), from the octet to nonet. In addition, new operators may appear too, as the unitary matrix 
u(x) is no longer traceless after the inclusion of the singlet. 

We recall the procedure to construct the U (3) xPT Lagrangian 18] before illustrating how to build 
the new operators in the resonance chiral Lagrangian. The key ingredient introduced to construct the 
U (3) xPT operators is the 1 /Nc expansion, in addition to the conventional expansion of momentum and 
the quark mass, which is usually named as the triple expansion scheme, i.e. 5 ~ p 2 ~ m q ~ 1/Nc- It 
also turns out to be useful when building resonance Lagrangian with the singlet 771 as an explicit degree 
of freedom 0TJ . 

As already mentioned previously, in this article we focus on the radiative decay processes with the 
types of VPj and P77, belonging to the odd intrinsic parity process describing two vector subjects 
(photon or vector resonance) and one pseudoscalar. Guided by the triple expansion scheme, the lowest 
order Lagrangian in the odd intrinsic parity sector is in fact the chiral anomaly formulated in the Wess- 
Zumino-Witten (WZW) action |42tl43j. with the order of 0(p 4 , Nc)- The relevant piece in our discussion 
can be written in the following way 

CWZW = -^f^ P a(^v"d^) , (5) 

where to get the photon field one needs to take = —eQ and the electric charge matrix of the light 
quarks with three flavours is Q = Diag{|, — |, — |}. 

The higher order Lagrangians can be categorized into two types: the higher order chiral anomaly 
pseudo-Goldstone Lagrangian and the Lagrangian with vector resonances. It is known that the higher 
order operators in the pure Goldstone Lagrangian encode the information of heavier degrees of freedom 
that have been integrated out. To avoid the double counting in RxT, the LECs of the higher order 
Lagrangian in the pseudo-Goldstone sector is usually assumed to be completely saturated by the heavy 
resonance states and thus the higher order operators in the pure pseudo-Goldstone sector can be dismissed, 
which works at least pretty well up to the 0(p 4: ) level in the even intrinsic sector [T^j. It is pointed out in 
Ref. [H] to fulfil this procedure it is necessary to use the antisymmetric formalism to describe the vector 
resonances. Though analogous analysis has not been carried out in the odd intrinsic parity sector, the 
resonance saturation assumption is utilized to construct the resonance Lagrangian in Ref. |25j and we 
generalize the discussion by including the singlet pseudo-Goldstone within the triple expansion scheme. 
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If the operator is written in terms of u(x), u(x) £ U(3), it obeys the canonical large Nq counting rules: 
terms with a single trace are of order Nq while one additional trace reduces its order by unity of I/Nq- 
The factor of ln(detw) also leads to a suppression of 1/Nc El HO]- The interacting vertex involving 
resonances and pseudo-Goldstones has the general structure at leading order of Nc 

d ~ (R 1 R2...R jX {n \v)}, (6) 

where x n (v 9 ) denotes the chiral tensor that only incorporates pseudo-Goldstone bosons and the auxiliary 
fields with the chiral order 0(p n ). For the odd intrinsic parity sector, it has been shown two types 
of vector resonance operators are relevant: (Vx^Hf)) and (VVx® (<p)) in the case of u(x) £ SU(3). 
When the singlet pseudo-Goldstone is taken into account, i.e. u{x) 6 U{3), two new operators with 
the same chiral counting order within the triple expansion scheme show up: (Yx if)) l n (det u) and 
(VV) ln(detu). The complete Lagrangians are found to be 

C V JP= 7^7 e^ pa ({V^J^}V a u-) 



+ ^({Va^JfR) 

- ic s M v ^e^ pa (V^fX)Hdetu), (7) 

£vvp= di e^dV^WnVaU*) 

+ id 2 e pvpa ({V^,vnX-) 

+ d 3 e pup(7 ({V a V^,VP a }u a } 

+ d4e P v P A{V a V^,Ve a }u a ) 

- idsM^e^iV^VnH^tu), (8) 

where we introduce tildes to the objects involving the Goldstone nonet to distinguish the one with octet. 
Comparing with Ref. [25], the new operators are 

0^j P = -ic 8 M y j| ff ^ CT {F^/f)ln(det«), 

5 VVP = -id 5 M^^e pupa (V^V pCT )Haetu). (9) 

We point out the above operators are only complete for the case with one pseudoscalar field. In the case 
with more pseudoscalar states, Ref. [25] has been generalized in Ref. [45J to include all of the relevant 
resonance operators that can contribute to the 0(p e ) xPT LECs in the odd intrinsic parity Lagrangian. 
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The relevant Lagrangian in the even intrinsic parity sector, describing the vector resonance and photon 
transition vertex, is [2] 

4 = ^(WD, (io) 

and the kinetic term for the vector resonance in the antisymmetric formulation reads [H] 

1 M 2 

where the nonet of the vector resonances resemble the flavor structure of the pseudo-Goldstone mesons 



P -vV° + + I ■ (12) 



For the vector resonances u and <f>, we assume the ideal mixing throughout this paper: 

2 /l 



The relevant R%T Lagrangian to our discussion can be summarized as follow 

C = C\VZW + £fcm + + £vjp + C-vvp ■ (14) 

4 Theoretical calculation of the radiative decay amplitudes 

In accord with the Lorentz symmetry, the general amplitude for the radiative decay V(q) — > P{p)i*{k) 
can be written as: 

iMv^Pj* = i e e^ pa e»e»* q p k a F v ^ Pr (Q 2 ) , (15) 

where ey and e 7 * denote the polarization vectors of the vector resonance and the off-shell photon respec- 
tively; the transferred momentum square is defined as Q 2 = —k 2 . In the case of the on-shell photon, one 
only needs to replace e 7 * with the on shell polarization vector e 7 and impose the real photon condition 
Q 2 = —k 2 = 0. Thus the decay widths of V — > P7 and P — > Vy are found to be: 

r(P -+ Vl ) = a{ M \ M fv f\F P ^ vr m 2 , (17) 

where a = e 2 /Air stands for the fine structure constant and the form factor Fp-^vy* (Q 2 ) can De defined 
in the same way as Eq. (|15p 

iMp(p)-+ V {q)-y*{k) = iee^ p(7 e l v e u 1 ,q p k CJ F P ^ Vr (Q 2 = -k 2 ) . (18) 
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The decay width of P — > 77 can be calculated 

T(P -> 77) = )va 2 Ml\F P ^ rr {0)\ 2 , (19) 

where the form factor Fp^.~ n *(Q 2 ) is defined in an analogous way as Eq. (|18p by replacing the vector 
resonance V with an on-shell photon. 

Next we use the resonance chiral Lagrangian described previously to calculate the decay widths and 
transition form factors defined above, which represent one of our main results in this work. We only 
consider the tree-level amplitudes here. The relevant Feynman diagrams to the radiative VPj* transitions 
are displayed in Fig. [TJ It is worth pointing out in the celebrated vector meson dominant (VMD) model 
only the type (b) diagram in Fig. [1] shows up. For the radiative transition P — > 77*, the Feynman 
diagrams are displayed in Fig. [2l 




Figure 1: Two types of Feynman diagrams for the processes VP^f* . 

Since the expressions for the various decay widths, such as V — > P7, P — > V7 and P — > 77, are 
rather lengthy, we relegate them in Appendix [A) The form factors of Fp^^*(s) and i 7 V_ > p 7 *(s) can be 
found in Appendix [Bj The kinematics and amplitudes of V — > -P7* — > Pl~l + , P — > Vj* — > Vl~l + and 
P — > 77* — > ^yl~l + , with the lepton I = e, /i are given in Appendix [Cj 

5 Short distance constraints from QCD 

Given R/\T being the dual theory of QCD in the resonance region, the couplings appearing in the 
resonance chiral Lagrangian can not be completely free, since there is only one parameter in QCD: 
Aqcd and the heavy quark masses. To fix the resonance couplings in terms of them requires solving the 
nonperturbative dynamics of QCD from the first principle, which is exactly the reason that the chiral 
effective field theory arises. Nevertheless, in the last decades, a sufficient way to implement the short 
distance feature of QCD has been developed, that is to match the operator product expansion(OPE) 
of Green functions of QCD currents which are of order parameter of the chiral symmetry breaking to 
the same quantity calculated within the resonance chiral Lagrangian [21\ \22\ [23l [2^1 [25] . Through this 
procedure, one could constrain sufficiently the resonance couplings in certain cases. 
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Figure 2: Feynman diagrams for -P77* 

The three-point Green function with vector-vector-pseudoscalar QCD currents has been studied in 
different works. Ref. [21] is devoted to the analysis of P — > 77 by taking into the higher order Goldstone 
chiral Lagrangian. The focus of Ref. [23] is to estimate the resonance contributions to the 0(p e ) LECs 
of xPT by using the vector formalism to describe the vector resonances, while Ref. [35] exploits the 
antisymmetric tensor formalism to incorporate the vector resonances. In Ref. [25], the VVP Green 
function has been analyzed by using the antisymmetric tensor representation for the vector resonances 
and the phenomenology study focused on the radiative processes of u — > 717 and ir — > 77. Our discussion 
in this section is devoted to the generalized study of V VP Green function given in [25] by extending the 
content of pseudo-Goldstone bosons from the octet case to the nonet one, i.e. u(x) € SU(3) — > u(x) £ 

The VVP Green function is defined as 

J d 4 x J d 4 ye^ x+ ^{0\T[V^x)V, b (y)P c (0)]\0) = d ahc e^p a q^ VVP {p\ q 2 ,r 2 ) , (20) 

where the flavor indices are a,b,c = 0, 8; is defined as r M = — (p + q)^; the vector and pseudoscalar 
currents are given by 



V-(x) = (^^)(x), P«( x ) = $i l5 ^)(.r). (21) 



Taking c = 0, one can get the singlet pseudoscalar current and Ilyyp can be defined correspondingly. 
Although the evaluation of the singlet Ilyyp and the octet n^p in the intermediate energy region could 
lead to different results due to the Ua(X) anomaly, their asymptotic behaviors coincide in the chiral limit 
and at leading order of 1/iVc 



(#) p 2 + g 2 + r 2 1 



lim n^p[(Ap) 2 , (Xq) 2 , (Xp + Xq) 2 } = lim U^ p [(Xp) 2 , (Xq) 2 , (Xp + Xq) 

A— >oo A— »oo 



[l + ( as )] + ( ). (22) 



2A 4 p 2 q 2 r 2 L ~\^n ~ K X 6 '' 
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Next let us focus on the calculation of the VVP Green function within RxT. Since the new operators 
given in Eq.© only contribute to the singlet Hyy P , the matching of the octet Green function ITyyp does 
not lead to any new constraints on the resonance couplings, comparing with the results in [25]. So we 
concei 
to be 



concentrate on the evaluation of the singlet VVP Green function Uyy P in RxT and the result is found 



n (0) (n 2 2 2, W)o [ AF 2 {d x - d z )r 2 + d 3 {p 2 + q 



+2V2 



F 2 \ v (M 2 - p 2 )(M 2 - q 2 )(M 2 - r 2 ) 
F v (ci + c 2 - c 5 )r 2 + (c 2 + c 5 - ci - 2c 6 )p 2 + (ci - c 2 + c 5 )q 2 



'My (M 2 - r 2 ){M 2 - p 2 ) 

+2V2 



F v (ci + c 2 - c 5 )r 2 + (c 2 + c 5 - ci - 2c 6 )q 2 + (ci - c 2 + c 5 )p 2 
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My (M 2 -r 2 )(M 2 -q 

32F v d 2 16V2F v c 3/ 1 1 ^ N c 



(My -p 2 ){M 2 -q 2 ) My { M 2 -p 2 M 2 -q 2 ' 8vr 2 (M 2 
+4\/3c 8 MyFy[— 3 J—, -+ ' 



-(M 2 -q 2 )(M 2 -r 2 ) (M 2 - p 2 )(M 2 - r 2 ) 1 

_ 9 /- d 5 F 2 M 2 \ 

(My - p 2 )(M 2 - ^(M 2 - r 2 ) / ' [Z6) 

We stress that the above expression for Ilyyp is worked out in the chiral limit. Due to the J7a(1) 
anomaly, the singlet pseudoscalar rji gains the non- vanishing mass Mq even at the chiral limit. Although 
the 771 mass Mq is suppressed by 1/Nc, its value is not a small quantity, which is even higher than the 
lowest vector resonance mass My [26]. Hence to take this effect into account we have introduced the 



non- vanishing mass for r\\ in the calculation of the Ilyyp function within RxT. 

Matching the result evaluated from OPE of the singlet VVP Green function displayed in Eq. (|22|) to 
the same quantity evaluated within RxT, which is given in Eq. (|23f) . leads to the following constraints 



4c 3 + ci = 0, (24) 

ci-c 2 + c 5 = 0, (25) 

„ „ N C My 

C 5" c 6 = CA _o 7^ . ( 26 ) 



64tt 2 v^F, 



V 



F 



2 



d 1 +8d 2 -d 3 = — T , (27) 



V 



N c M 2 F 2 V3M V „ V2M 2 . 

+ ^2" ^ C « ~ U 77 c l • ( 28 ) 



647T 2 F 2 8F 2 Fy 8 MyFy 

We find the above constraints are all consistent with the ones given in [25], except for the relation of d 3 . 
The consistent condition for the results of d 3 from the octet and singlet cases requires 

h= -^m h - (29) 

So in the numerical discussion, we will take this constraint for c%. 

There is another well known approach to address the features of form factors corresponding to exclusive 
processes of QCD with high momenta transfer, which was developed within the parton description scheme 
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for the hadrons in Ref. The relevant one to our present discussion is the photon meson transition 
form factor F 7T ^(Q 2 ), which can match the form factor defined in Eq. (|15p by applying the time reversal and 
replacing the vector resonance with the on-shell photon. Although different approaches have predicted 
different asymptotic behaviors of F^^Q 2 ) at the order of Q~ 2 , for example 



F nj (Q 2 -> oo) 

F-K-yiQ 2 ->■ OO) 



F 

Q 2 ' 
2F 



3Q 
F 



2 



F^Q^oo) = __ ; (30) 

which are noted in Refs [46j[47|[48j respectively, they agree at the order of Q°, indicating the form factor 
behaving smoothly at the high momentum transfer, i.e. F n ^(Q 2 — > oo) — > 0. So the most conservative 
constraint would be just to impose the vanishing condition for the constant term in the form factor, i.e. 
to demand the coefficient of Q° being zero. The explicit expression for the form factor of ^77* can be 
found in Appendix [C] and the corresponding high energy constraint from the order of Q° is 

c\ - c 2 + c 5 = 0, 

~ ~ JV T , N C My 

C5 " C6 = 7m 3 + ^fr^v' (31) 

and the corresponding constraint from the form factor of ojtt^* , which is also given in Appendix [Cl leads 
to 

c 1 -c 2 + c 5 = , 

c 5 -c 6 = -^4- (32) 

Combining Eq. (|3ip and Eq. ()32p . we have the following relations 

ci - c 2 + c 5 = , 

N C My 



C5 - C 6 



64vr 2 



v 



N c M v 



d3 = -6^4 • (33) 

When discussing the high energy constraints, we take the chiral limit and U(3) symmetry for the vector 
resonances. In this case, the physical states r\ and r/ are represented by the flavour eigenstates % and 
771, respectively. So in the chiral limit, the constraints from the other processes with octet pseudoscalar 
mesons lead to the same results, while the process involving the singlet pseudoscalar rj\ could lead to 
different results as we keep the Ua(1) anomaly effect, i.e. the non-vanishing r\\ mass, in the deriving the 
high energy constraint. The explicit result from the analysis of the 7/177* form factor is 

c 1 -c 2 + c 5 = 0, 

N C My 



C5 - c 6 



64vr 2 ^2F, 



v 



N c M v V3M V „ V2M$ 



647T 2 F 2 Fy 8 MyFy 



c 8 - j-r-^ci , (34) 
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To demand the consistency of the results from ^77* and 77177*, we arrive at the same constraint in 
Eq.p}. 

Notice that the only inconsistency from the OPE and form factors is the relation of d% , although fewer 
constraints are obtained in the analysis of form factors. This observation has been confirmed in other 
processes |19t [20] . The fact that only one multiplet of resonances was unable to fulfill all the high energy 
constraints from OPE and the corresponding form factors was already noticed [24 ^ I23[ 139]. In this work, 
we will take the result for d% in Eq, ([33p obtained from the high energy constraint of form factor, which 
is more related to the processes we are discussing. To reduce the free parameters as many as possible in 
the phenomenology discussion, we also exploit the constraints in Eq. (|24p and Eq. ()27p obtained from the 
OPE analysis. In summary, the high energy constraints we are going to use in the following discussion 
are those in Eq.fl23J), Eq.flZZD, Eq. ljgS]) and Eq.fl33|). 

6 Phenomenology discussion 

Although we can fix some parameters through the short distance constraints from QCD in the previous 
section, some of the resonance couplings appearing in the decay widths given in Appendix are still 
unconstrained. So we need to fit the unknown resonance couplings, such as 03,^2 and d$, together 
with the Tj — rf mixing parameters 9q,8s,Fq and F%. For the mass of vector resonances in the chiral 
limit My, one can safely estimate its value by M p , the mass of /o(770) [17] , While for the parameter 
Fy, describing the transition strength of the neutral vector resonances and photon, its value is still 
a somewhat controversial subject and several solutions have been proposed in different works. In the 
pioneer work to discuss the high energy constraints in RxT [33]. Fy = V2F was predicted by combining 
the high energy constraints from the pion vector and axial-vector form factors within the minimal RxT 
Lagrangian at leading order of 1/Nc, while in the next-to- leading order analysis of the pion vector form 
factor Fy = \^3F is updated in Ref. [50J , which has also been confirmed in the study of radiative tau 
decays [20] and the partial wave analysis of tttt scattering [51} 152] . Phenomenology determinations of 
Fy = 147 MeV and 180 MeV have been used in r — > VPv T [18] and r — > KK-k [19] decays respectively. 
By estimating the pseudo-Goldstone decay constant at chiral limit by the pion decay constant F n , one 
has 

F v = V2F = 131 MeV , Fy = V3F = 160 MeV . (35) 

Apparently, a more precise value for Fy is needed in our discussion, since the physical processes we are 
discussing are mainly the radiative decays of vector resonances and precisely Fy describes the interaction 
of vector resonances and photon. Thus we decide to free Fy and fit its value in our program, in such a way 
we could predict a more reliable value, as Fy is rather sensitive to the processes we are considering. To 
exploit the high energy constraint for eg in Eq. (|29p . the value for the C/a(1) anomaly mass Mq is needed, 
which has been reviewed in Ref. [2B] and also recently determined in Ref. [52]. We use the average value 
Mq = 900 MeV from the two mentioned references throughout. For the other unmentioned inputs, unless 
an explicit statement is given, we will take the corresponding values from Ref. [12] . 

As advertised previously, although the main contribution of our current work is to determine the rj—rj' 
mixing parameters in a more reliable theoretical framework, we include the relevant radiative processes 
without rj and rf into our discussion as well. By performing the \ 2 fit, we can determine the unknown 



11 



resonance couplings and rj — rf mixing parameters: 



F 8 = (1.37 ± 0.07)^ , F = (1.19 ± 0.18)^ , 

(9 8 = (-21.1 ±6.0)° , O = (-2.5 ±8.2)°, 
Fy = (136.6 ± 3.5)MeV , c 3 = 0.011 ± 0.016 , 

d 2 = 0.086 ± 0.085 , d 5 = 0.36 ± 0.40 , (36) 

with x 2 /d.o.f = 64.0/(70 — 8) = 1.03. For the various decay widths, we summarize the experiment data 
and the results from our fitting program in Table [1] for the processes without r\ and rj' and in Table [2] 
for those involving r/ or rf. To visualize the results, we plot the numbers in Tables [T] and [2] in Fig. [3] 
where one should notice we have scaled different decay widths to a proper range in order to show them 
in one figure. The resulting plots for the form factors of t?77*, 7/77* and (pr]j* are given in Figs. HI 
[5] and El respectively. The error bands shown in the plots and the errors of the parameters in Eq. (|36p 
correspond to the statistical uncertainties at 2 standard deviations [53]: n a = (x 2 — XoVV^Xoj with Xo 
the minimum x 2 obtained in the fit and n CT the number of standard deviations. 





Exp 


Fit 


Theo (F v = 160 MeV) 


Theo (F v = 180 MeV) 




757 ± 28 


731 ± 37 


533 


421 


r p n^ 7r o 7 


89.6 ± 12.6 


76.0 ±38 


55.4 


43.8 




116 ±12 


113 ±6 


83 


65 


^uj— >ire~~ e+ 


6.54 ±0.83 


6.64 ±0.33 


4.84 


3.83 




0.82 ±0.21 


0.66 ±0.03 


0.48 


0.38 



Table 1: Experimental and theoretical values of the various decay widths without r\ and 77'. The experi- 
ment data are taken from |12| . All of the values are given in units of KeV unless specified. The results 
from our fit are listed in the column Fit and the error bands of the widths are calculated by using the 
same parameter configurations that we use to get the error bands for the parameters in Eq. ()36p . To 
show the relevance of Fy, we have given another two theoretical predictions for the various decay widths 
in the last two columns by taking Fy = 160, 180 MeV. 





Exp 


Fit 


p 

x UJ — ¥7]~f 


3.91 ±0.38 


5.05 ±0.36 


TpO — ^777 


44.8 ±3.5 


41.6 ±3.2 


»r/7 


55.6 ± 1.6 


55.3 ±2.5 




0.265 ±0.012 


0.270 ±0.021 


p , 

A rj'—^wy 


6.2 ± 1.1 


7.4 ± 1.0 


A 77—^-77 


0.510 ±0.026 


0.481 ±0.038 


x 77'— >77 


4.30 ±0.15 


4.25 ±0.21 




(8.8 ± 1.6) x 10- 3 


(8.0 ±0.6) x IO- 3 


^77— >jfi~ /J.+ 


(0.40 ±0.08) x 10~ 3 


(0.38 ±0.03) x 10~ 3 


p 

1 7/'— >7/J, — /J,+ 


(2.1 ±0.7) x 10~ 2 


(1.8 ±0.1) x 10~ 2 


p 

1 0— )-77e~ e+ 


0.490 ±0.048 


0.464 ±0.021 



Table 2: Experimental and theoretical values of the various decay widths involving rj and rj' . The 
experiment data are taken from [12] . All of the values are given in units of KeV unless specified. The 
error bands of the widths are calculated by using the same parameter configurations as used in Table [1] 
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Different modes 



Figure 3: Different decay widths. We have scaled different decay widths into a common region in order 
to show them in one figure. For the values before scaling, see the numbers in Tables [1] and The open 
squares denote the central values with the best fit given in Eq. (|36|) . and the shaded area correspond to 
the error bands generated by the parameter configurations explained in the text after Eq. (|36p . Note that 
using the high energy constraints we can completely predict the decay widths of ir — > 77 and tt — > je~ e + , 
which are in good agreement with the experimental data, and here we include these two processes just 
for completeness. We do not include rj' — > p-y in the fit, since PDG [12J also includes the background part 
from rj — > ixix^f to determine the width for rj — )■ pj. Nevertheless due to the dominant decay channel of 
the p resonance is irir, our prediction for 7/ — > pj agrees with the one from PDG [12J. 

Several remarks about the fitting results are in order. We comment them as follows. 

1. The first lesson we can learn from the results in Eq. (|36p is that 63, d<i and cis carry huge error 
bars. Nevertheless we find there exist strong correlations among these parameters. We plot the 
correlations of ^2-^5, c^-di and 63-^5 respectively in Figs.©, ([8]) and ([9]), where the same parameter 
configurations have been used as we exploit in evaluating the error bands in Eq. (|36p and Figs.([3]l6|). 
A very strong linear correlation between d.2 and is observed, as one can see in Fig.©. While the 
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Figure 4: The form factors of rj — > 7*7. The solid line (red) denotes the result from the best fit and the 
shaded area correspond to the error bands. Sources of the different experiment data are: solid squares 
[MIES]) open squares [56], open circles [57], solid triangles [58], open triangles [59]. The separated 
figure in the upper part is the close-up of the main plot in the region of s > 0. 

correlations of 03-^2 and £3-^5, as shown in Figs. ([8]) and ([9]), are not as strong as ^2-^5. A correlation 
between C3 and cfo (in fact they are C3 and di from the pseudo-Goldstone octet Lagrangian) has 
been revealed in a preliminary analysis of r — > r K r KT)v T decay [BT], where r\ particle is treated as 
the pure octet rj$. The parameter space found in the previous reference covers most of the space 
shown in Fig. ([5]), but the correlation relation in jfJJJ is with opposite sign of the relation we find 
in the present work. This indicates that the combining study of the radiative decay processes of i] 
or rl and the r decays involving rj or 7/ may help us pin down the resonance parameters, such as 
£3 and c?2, which deserves a future work. For the remaining parameters in Eq. (|36p . we do not find 
significant correlations among them. 

2. The value of Fy. As one can see in Eq. (|36j) . our analysis favors a smaller value for Fy, comparing 
with the values used in [T2] and |19j . By using the high energy constraints mentioned in previous 
section, the processes appearing in Table [T] are solely determined by Fy. Thus those radiative 
decay processes without r\ and rj provide a tight constraint on the value of Fy. To clearly show 
the relevance of Fy in those channels, we give another two predictions by taking Fy = 160 and 180 
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Figure 5: The form factors of rf — > 7*7. The solid line (red) denotes the result from the best fit and the 
shaded area correspond to the error bands. Sources of the different experiment data are: solid squares 
[Ml EH], open squares [56] . open triangles [59], solid circles [60J. 

MeV in the last two columns in Table [TJ We also use these 5 processes in Table [1] to perform a fit 
to get the value of Fy, and the result is Fy = (134.05 ± 2.17)MeV, with x 2 /d-o.f = £z? = 0.36. 
This result for Fy is in perfect agreement with the global fit in Eq. (|36p , If one takes the value of 
Fy = 134.05MeV and fits the other 7 unknown parameters by using the rest of experimental data 
as we used to get Eq. ([3"6"j) , the fit results turn out to be quite similar to the ones we show in Eq. ([33]) , 
as expected. 

3. The rj — rj' mixing parameters. For Fg, it can be completely fixed by the ratio of Fx/F n in the 
next-to-next-to-leading order within the triple expansion of large Nq xPT, which yields a rather 
reliable prediction F$ = 1.34^ [11] . However at the same order, there exist several unknown 
LECs for the predictions of Fq and #0 — 9s, which prevents the precise predictions for their values. 
As one can see in Eq. (|36p . our result for F$ agrees with the xPT prediction. In Ref. [10] . Fq was 
determined in the process P — > 77 at next-to-leading order by ignoring the chiral symmetry breaking 
operators. Assuming OZI violating coupling in the next-to-leading order to vanish, Fq = 1.25F n can 
be derived [1Q|, 131] , Fq can be also evaluated in the standard way to include the chiral corrections 
from loops and LECs in the calculation of the axial- vector current matrix element [11] , In this 
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Figure 6: The form factors of <f> — >■ 777* |57j . The solid line (red) denotes the result from the best fit and 
the shaded area correspond to the error bands. 

case, two additional OZI violating couplings appear and if one assumes those couplings to vanish, 
Fo ^ -Fjt can be predicted. So a reliable determination of Fq could help us better understand the 
somewhat inconsistent results from the two approaches. As an improvement, we have included not 
only the OZI violating operators but also the chiral symmetry breaking ones in R/\T to determine 
the value of Fo in P — > 77 processes. Our result shows the inclusion of the OZI suppressed and 
higher chiral symmetry breaking operators in the calculation of P — > 77 does not change the result 
of Fq = 1.25F n very much, as obtained in |1U] , This also indicates that the ignorance of the OZI 
suppressed operators in the calculation of the axial-vector current matrix element, which leads to 
F ~ F w , is not a good approximation. About the mixing angles, our result 9$, = (—21.1 ± 6.0)° 
agrees well with the results in literature, see Table 1 of Ref. [26J. While for the mixing angle 8q, 
our analysis reveals that a huge error accompanies this parameter, as shown in Eq. (|36p . This may 
be viewed as a source that why rather different results have been obtained for #0 |26j. 

By using simple parameterizations of the VP7 and P77 vertexes that only consist of constant terms 
(independent of the quark masses and momenta), the rj — rf mixing parameters have been explored 
in various V — > P7 and P — > 77 processes and one can sec Ref. [26j for a comprehensive analysis 
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Figure 7: The correlations between di and d$. The solid line (black) corresponds to d$ = 4.4^2 — 0.06. 

of different results. Our current work confirms the validity of the two- mixing- angle description 
in a more general framework for VP7 and P77 interaction vertexes. Within the current theoret- 
ical framework, the various data from different radiative decay processes can be simultaneously 
incorporated and our x 2 is clearly better than that in Ref. |38| . 

4. An additional interest of our work is to test the resonance saturation assumption for the couplings 
relevant to the chiral anomaly induced processes P —> 77, with P = ir,r],r]'. The next-to-leading 
order odd intrinsic parity operators in the chiral Lagrangian with only pseudo Goldstone bosons 
can be categorized into two parts: 0(N^p e ) and 0(N^p 4 ) \21\ 162] : 

L ( l = & W<*-/f /+*> ~ hs^(V x f^{ff, u"}) + h ) V60o. (37) 

The higher order LECs appearing in the pseudo-Goldstone chiral Lagrangian encode the high energy 
dynamics of the underlying theory. In the resonance saturation approach, it is assumed that the high 
order LECs are completely saturated by the resonances and thus one does not need to include extra 
pure higher order pseudo-Goldstone operators in the resonance Lagrangian. This assumption has 
been proven to be very successful for the 0(p 4 ) xPT LECs in the even intrinsic parity sector [14] . 
In the present discussion, we shall check this assumption in the odd intrinsic parity sector. By 
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Figure 8: The correlations between £3 and d%. The solid line (black) corresponds to d% = 5.6C3 + 0.06. 



integrating out the vector resonances in the resonance chiral Lagrangian introduced in Sect. [3j we 
have the following predictions 



'2 



Fv F 2 
1 (ci + c 2 + 8c 3 - c 5 ) + {di + 8d 2 -d 3 ), 



V 



F V ,- , Fy , 

(C5 - C 6 ) + 7TT7r«3 



v^M 3 







6\/2M, 3 



2My 



(ci + c 2 - c 5 ) 



F y c 8 , F 2 M 2 



2 v / 3Af 1 



+ 



\2My 



[d x - d 3 ) + 



Fyd 5 



(38) 
(39) 
(40) 



where , ^ coincide with the results in Ref. [25J and the result for ^3 is a new result to our knowl- 
edge. Taking into account the high energy constraints of the resonance couplings in Eqs. (|24p (|27| )(|29p 
and (f33|) . we have the simplified predictions 

F 2 zv _ N C L y_F 2 M 2 IFyMl - | F 2 ^ 



64M y 



64vr 2 M 2 ' 



F 2 Mq 2F 2 M 2 

«2 + 



96M y 



3M y 



2V6M 2 



(41) 



The fitted results in Eq. ([36|) yield 

l x = 0.37 x 10~ 3 GeV" 2 , t 2 = -8.0 x 10~ 3 GeV" 



(0.83 ± 3.02) x 10" 



(42) 
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Figure 9: The correlations between £3 and d§. The solid line (black) corresponds to d§ = 25c3 + 0.2. 



where the error of £3 is estimated by taking the errors from the fit results in Eq. (j36p . 

Now we can calculate the form factors of P — > 77 using the pseudo-Goldstone Lagrangian in Eq.([5]) 
and Eq. ([37j) . For tt — > 77, the k% operator is irrelevant and our result is the same as the one in 
Ref. [25] 

zp _ pWZW _i_ rpti /,iq\ 

7r vyy — 7r — ^'yy ' r Tv^jj j \^°) 

where F^^^ denotes the contribution from the WZW Lagrangian in Eq.([5]) 

i?wzw = l _ (u) 

r ^77 4vr 2 F^ ' K ' 

and F%L, denotes the contribution from the t\ term 

64 



FlU^ = ^rrnlh . (45) 



Note that the £2 operator does not contribute to the considered process. For 77 — > 77, we have 

p — pWZW 1 pti , pk 3 (AP.) 
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where 



F. 



wzw 



1 . cos 0q sin 8s 



cos (0 O - 8 ) 4V3vr 2 v F 8 



), 



and F*!^, F^s denote the contributions from iii^3 operators respectively 



COS (fo 



r 64 2 „ 2 x cos #0 . 

■\ — F 7m. -Mif — h ^- 



128^2\ 2 ^ 2 -sin^s ~ 
(2771^ + m#) — 



F 



(47) 



(48) 



Similar result for rf 



where 



- 1 T]— >77 



1 64\/6sm0 8 



cos (0 O - ^8 ) 3F 
77 is found to be 

p , _ pWZW , F k , pfc 3 



"rf— >77 



cos (0 O - ^8 ) 4v/3tt 2 F f 



1 sin 0q _|_ cos 08 



, n r 64 2 , 2 ^sin0 o 128\/2\ 2 , 2^ os ^r 
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(49) 

(50) 
(51) 

(52) 



rf — >77 



64^/6 008 8 



cos (yo 



3F 



(53) 



In Table [3l we show different contributions to the form factor Fp_s. 77 using the predictions for t\ 
and k% in Eq. (|42p . For the tt — > 77 process, the chiral symmetry breaking effect is rather tiny, about 
1%, since the chiral correction is proportional to m^/My. Hence the leading order contribution 



from the WZW Lagrangian overwhelmingly dominates the decay width of it — > 77 [25, 36J. For 
57 — > 77 and rf —> 77, the WZW term can also give rather close results to the experimental values, 
as one can see in Table [3] the next-to-leading order correction is only at most 14% of the WZW 
term. So our current calculations confirm the validity of the triple expansion scheme for the odd 
intrinsic parity pseudo-Goldstone Lagrangian, i.e. the WZW contribution plays the dominant role 
in the P — > 77 processes. We point out that this conclusion is based on the fact we have used the 
fitted results for the mixing parameters in Eq. (|36p . 



Table 3: 
in Eq. 





rpWZW 


77*1 


pk 3 


^WZW+h+ks 
^P^-77 


r P-S>77 


7T — > 77 


-0.274 


0.002 





-0.272 


-0.275 ± 0.070 


rj — > 77 


-0.265 


-0.007 


0.015 


-0.256 


-0.272 ± 0.070 


7/ — > 77 


-0.365 


0.011 


0.039 


-0.315 


-0.342 ±0.012 



The predictions of different contributions to the form factors Fp_> 77 using the parameter values 
All of the values are given in units of GeV -1 . 
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5. We can predict the decay widths of p — > ire + e , rf — > je + e and <p — > rjp + p by using the results 
from the global fit Ea.flMl): 

Tp^e+e- = (3.42 ± 0.17) x lO^KeV, T v ,^ e+e - = (8.85 ± 0.48) x 10" 2 KeV 
r^ w+At - = (2.22 ± 0.13) x 10- 2 KeV, (54) 

These predictions are below the upper limits given in PDG |12| . Hence our results can provide a 
theoretical hint to the future experimental analysis on these channels. 

7 Conclusion 

In this work, we complete the resonance chiral Lagrangian in the odd intrinsic parity sector ( VVP type) 
by including the singlet rj\ field as the dynamical degree of freedom. We exploit this Lagrangian to study 
radiative decay processes: P — > V7, V — > Pj, P — > 77, P —¥ V — > Pl + l~, as well as the form 

factors of r] — > 77*, rf — > 77*, <fi —> 777*. The two-mixing- angle scheme is used to describe the rj — rf 
system in the discussion. By imposing the proper short distance behavior of QCD, we can fix several 
combinations of the unknown resonance couplings. The remaining free resonance parameters, together 
with the mixing parameters Fg , Fq , #g , #0 are determined through fitting the various experimental data. 
We have shown the resonance chiral Lagrangian can provide a systematic theoretical framework to handle 
the various radiative decay processes involving the resonance states and simultaneously accommodate the 
various experimental data. Thus we believe the rj—rj' mixing parameters resulted from this analysis should 
be rather reliable. 

By integrating out the resonance states in R%T, we predict the higher order low energy constants in the 
odd intrinsic parity pseudo-Goldstone Lagrangian. We conclude the WZW contribution from the leading 
order dominates the processes of P — )■ 77, with P = Tr,r] and 77', by using the 77 — r( mixing parameters 
from our current analysis. We have also predicted the decay widths of p — > ire + e~ , rj' —> je + e~ and 
<p — > r]p + n~ , which may shed light on the future measurement for these three channels. 
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A The decay widths of V — > P7, P — > V7 and P — > 77 

The various decay widths from different processes are given below 

1 Ml -ml 



T(o; — > 7T7) = -a( 



^3 



3 v 2M L 
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( 2/2 

+ Fn M^ M 2 + 8d 2 - d,)ml + d 3 M*] }' , (55) 
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V3M y M (~FT + F ^ [ ' ( ° 2 " Cl + C5 " 2cg ) + m ''( C2 + Cl " C5 ) + 8c3m '] 

4F y , sin 0p \/2cos08 U ~ 2 2 ~ 2 ~ 2] 
- ^ M3 (-^r- + ) L d 3(M w - ny) + im,, + 8d 2 m w J 

,COS0 8u 4x/2My_ 8FyMy j^] 2 



r(7r -)• 77) = ^na 2 Ml 

f iVc 4ViFV / 1 , l u -^ xr ^ 2 

1 " 12^ " sfJ^^m] + m^ {ci + C2 + 8cs - C5)m - 

8F 2 , 1 ^ 2 



V (172172 Kdi + fcfc-^K !> - (63) 



IF^MIMI' 
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1 2 3 1 f N c ,cos6 2\/2sin6» 8 , 



W f " 4 "cos (0 o -0 8 ) I 12^vr 21 F 8 F y 

^/2sin 6 8 s 

~vm; { ^i + ^ [m " (C2 + C1 " C5) + 8csmJ ( ^ ^ } 

\QF V 1 r 2 ~ \ i q~ /'o 2 aM/^cos^o . sin0 8 . 



7T { Mj + 9M* ] ^ ~ ^ + M2m ^ { — s 

8y/2F$ 1 r~ 2 t 2 , fi t ,„ 2 2M/>/2cos0o . sin0 8 , 
9^3 ^ " ~ 17 8d2 ( 2niK ~ ( ^ + 



+ (-^)[- 8 y2M,^ 8 (^ + ^ + ^: 



+ 8F^-i 5 (-^ + ^ + s |j)]} 2 . (64) 



1_ 2 3 1 f Afc /Sin^o 2\/2cosw 8 s 



W 7 4 "cos (0 o -0 8 ) I 12^7T 2V F 8 F ' 

4\/2Fy 1 1 u 2 /~ ^~ ~\,e~ 21/ sin ^0 , V2COS08, 

"7i^ ( mj + ^ } K ' (C2 + C1 " C5) + SC3m * ] + — ) 

16 F v 1 r 2 ,~ . ~ ~ x , c ~ fo 2 2M / ^singp cos 6> 8 > 

+ K ' (C2 + C1 " C5) + 8cs( mx " mJ] ( ^ ^ } 

, AF V ( 1 1 \ r T 2 7 2 i qJ 2U sin #0 . \/2~COS(9 8 , 

+ V5 ( mJ + [iim " ~ + 8 2 '"- ] + i^i 



+ 8^M?*(- j | + ^ + s |j)]} 2 . (65) 
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B The form factors of — > 777*, 77 — > 77* and rf — > 77* 

The definition of the form factor is given in Eq. (|15p and the explicit forms for different processes are 
reported below 



<fi—>vy 



1 



COS (Oo — U£ 

Ay/2 ,y/2 cos 6*o sin^gx 

+ — =— ) x 



x [M|(c 2 - h + c 5 - 2c 6 ) + m 2 v (c 2 + h - £5) + 8c 3 (2m^ - m^) + (ci - c 2 + c 5 )s] 

8i?l/ n / wV2cos0o . sin0 8u ~ 2 2 - 2 ~ 2 2 - 

( h — ^— J ^(M^, -m v + s) + dim^ + 8d 2 (2m K - mj 



, sin 5 8 u 8My . 8y2/ r .Uf , 
+ ( " IT) [ - 3Ml C8 + 3M, <* 5 ] f ' (66) 



where the definition of Dr(s) is 



M R - s- iM R T R (s) 



For the narrow-width resonances u, (f), we use the constant widths for them in the numerical discussion. 
For p resonance, the energy dependent width is constructed in the way introduced in [63j : 

I» = 9^§2-[^(* " 4m*) + ^0( S - Am 2 K )), (68) 



where op = J\ — Amp/s and is the step function. 
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The form factors for 77 — > 77* and rj' — > 77* are 

F ^ As) = cos(el-e 8 ){ 

N c ,cos<9 2A/2sin6» 8 , 



12^3^ 1 F 8 " ~ F 
2V2FV 1 J_ . 



r 2/~ , - - x , o~ 2 , /- - , - x IfCOsOo y/2 sin fl 8 ^ 

[m v (c 2 +ci - c 5 ) + Scsm^ + (ci - c 2 + c 5 )sj (— — J 

■^[ D ,(.) + ^W]x 

r 2 /~ , ~ ~x , o~ 2 , /- , ~ - ~ n 1/ cos O sin 6> 8 , 
m (c 2 + ci - c 5 ) + 8037^ + (c 2 + c 5 - ci - 2c 6 )s (— ) 

^8 ^0 



8Fy , 1 , 
H 7= ( — o) x 

r 2 2 2 l / cos #0 sin 08 . 
m (C2 + ci - c 5 ) + 8c 3 (2m K - mj + (ci - c 2 + c 5 )s If h — ) 

^8 ^0 

8F v n / w \/2cos0 o , sin<9 8 , 
[m^(c 2 + ci - c 5 ) + 8c 3 (2m|- - m^) + (c 2 + c 5 - ci - 2c 6 )s ] 

^W] M + d 3 (s - ml) + M 2 {2m\ - + ^) 

+ ( - ^) [ - 4V2M vFv c 8 (D p{ s) + ^ + l Dw(a) + _L + |^ (a) + _|_) 
+8F y M y d 5 (^^(,) + ^(.) + D^)^)] } , 
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cos (6»o - 6*8 ) 



N c / sin 6»o 2\/2cos6» 8 s 



_V2^V J_ J_ 

r 2 2 n /■ srn $0 COS #8 \ 
|m T; ,(c 2 + ci - c 5 ) + Sc^m^ + (ci - c 2 + c 5 )sj (— — H — j 

r 2 2 n / sm ^0 COS #8 \ 

m ,(c 2 + ci - c 5 ) + 803771,,. + (c 2 + c 5 - ci - 2c 6 )s f — 1 ) 

8F V . 1 w >/2 sin #0 cos $8 \ 

[771,2, (c 2 + g x _ g 5 ) + 8c 3 (2m|- - m 2 ) + (ci - c 2 + c 5 )s] 
8F v n /■ w V2sin0 o cos6»8x 

+ ^Uh d * {s){ — «T )X 

[m^/(c 2 + ci - c 5 ) + 8c 3 (2m 2 K - ml) + (c 2 + c 5 - Ci - 2c 6 )s ] 
4Fy r 1 _ . , 1 / \i r 7 2 7/ 9\ o7 o-i/Sinflo \/2cos^8\ 

[^W] + - *%) + 8d 2 (2< - tti 2 )] - ^) 

+ 8^M^ 5 (^^ p ( S ) + ^ + ^W^)] } • (70) 

C The decay widths of V -> P/"Z+, P /+ and P -> 7/ /+ 

The kinematic variables used in this sector are defined below. 
• For V(p) -> l-{k 1 )l + (k 2 )P{q) 

k = ki + fc 2 , 
(fci+/c 2 ) 2 = s, 
(g + fea) 2 = t. 

(71) 

The general decay amplitude is 

T(V -> rZ+P) = e Vp^*V ufaJVyk) ^- Z+P (s, i). (72) 
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To calculate the decay width, one needs 

2mf s + 2m 2 [mp + M v - M v s - m 2 P {2M v + s) - 2s t] 

+s [mp + M v + s 2 + 2s t + 2t 2 - 2m 2 P (s + t) - 2M v (s + t)] . 

Then the decay width is found to be 

I l l r(M v -mp) 



V ^ Pl l+ 3 8tt 3 32 My J 4m 2 



1 AMv-vapY ri mM 

m- y L d - L * |T(M)|! 



where 



M 2 +m 2 p + 2m 2 -s ^ s{s - 4m 2 )[s - (My + m P ) 2 ][s - {My - m P f 

trriin — ~ 



2s 



Ml + m 2 +2m 2 -s J s ( s ~ 4m U s ~ ( M v + m P ) 2 ][s - (My - m P f 



' + 



For P(q) -»• l'(k 1 )l + {k 2 )V{p) 

k = k\ + k 2 , 

(h + k 2 ) 2 = s, 

( P + k 2 ) 2 = t, 



2s 



The general decay amplitude is 

T {P -> = e 2 e lxl/pa e^ P a u(k 1 )^v(k 2 )Ap^ l - l+ y(s,t). 

To calculate the decay width, one needs 

\e^ P ae;kfp- u(h)^v(k 2 )\ 2 = 

2m\s + 2m 2 [mf> + M v - M v s - m 2 P (2M v + s) - 2s t] 

+s [mp + My + s 2 + 2st + 2t 2 — 2m 2 P (s + 1) - 2M v (s + t)] . 

The decay width is 
where 

M 2 + m 2 + 2m? - a \J < s ~ 4m U s ~ ( M v + m P ) 2 ] [s - {My - m P ) 2 } 

2 2s 
M 2 +m 2 P + 2mj - s \/ s ( s ~ 4m U s ~ ( M v + m P ) 2 ][s - {My - m P ) 2 } 
2 + 2a 



tmax 
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For P{p)^l-{k 21 )l + {k 22 ) 1 {k l ) 



k 2 = k 2 i + k 22 , 
(k 21 +k 22 ) 2 = s, 
{k 1 +k 22 f = t. 



I W e i*i fc 2 u(k 21 )^v(k 2 2) | 2 = 

2mf s + 2m 2 (m P - m 2 P s -2st) + s[m% + s 2 + 2s t + 2t 2 - 2m|>(s + i)] . 



The decay width is 

cm?, ctmax 



I f m p rtmax 



where 



m p + 2m 2 — s 
tmin = — ■ 



(nip — s)ijs(s — Am 2 ) 



2 2s 



2 p + 2m 2 -s (m 2 p-s)^s(s-4m 2 ) 



m 

tmax = — — h 



2s 



The explicit expressions of the decay amplitudes are given below 

T(u -> ttI-1 + ) = -e 2 e Wff ^V u(h)^v(k 2 ) J | 

2^/2 



(81) 



The general decay amplitude is 

T(P -> r i+7) = e Vp^i*^ u(fc 2 ih'M*22) A P _ w -, +7 (a, t). (82) 
To calculate the decay width, one needs 



(83) 



(84) 



(85) 



F^MyM u 



[(ci + c 2 + 8c 3 - c 5 )m 2 + (c 2 + c 5 - ci - 2c 6 )M 2 + (ci - c 2 + c 5 )s] 



^ p(s) + M ~ 2 " d ~ 3)m * + d ~ 3(M * + s)] } ' (86) 



T (p 7rrZ +) = _ e 2 e^e^V ^1)7^2) ^ | 

2^/2 



3F T M V M. 



p 



[(ci + c 2 + 8c 3 - c 5 )m 2 + (c 2 + c 5 - ci - 2c 6 )M 2 + (ci - c 2 + c 5 )s] 



+ sSp DUs) + 8d ~ 2 " Js)m ' + Js(Mp2 + s)] } ' (87) 
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T{<j)^r]ri + ) 



= — e e 



e v v k p f u{k{]^v{k 2 ) - (a) 



(88) 



r(7r^ 7 r/ + ) 



N c 2V2F V , 1 1 , r ,„ _ o _ _ . 2 ,„ _ _ , , 
"12^ " 3F^^W p + M?) t (Cl + C2 + 8C3 " C5)m - + (C1 " C2 + C5)S ] 




(89) 



T{r t ^ 1 l-l + ) 



e 3 e fiupa eik^k% u(k 2 i)^v(k 2 2) -i^ 77 * ( s ) 



(90) 



pvpu^xKkl u(fc 2 i)7 M f (£22) -i^'-> 77 * (s). 



(91) 



where F^yj* and F^^^* are given in Eqs. ([69"j) and (fTOj) . 
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